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Abstract 



In this article, we calculate the scalar form-factor of the proton in the 
framework of the light-cone QCD sum rules approach with the three valence 
quark light-cone distribution amplitudes up to twist-6, and observe the scalar 
form- factor a(t = —Q 2 ) at intermediate and large momentum transfers Q 2 > 
£\\ ■ 2GeV 2 has significant contributions from the end-point (or soft) terms. The 

numerical values for the a{t = —Q 2 ) are compatible with the calculations 
from the chiral quark model and lattice QCD at the region Q 2 > 2GeV 2 . 

o ' 

PACS : 12.38.-t, 14.20.Dh, 13.40.Gp 

Key Words: Sigma term, Light-cone QCD sum rules, Scalar form-factor 

O 

"Sh! 1 Introduction 

The pion-nucleon sigma-term E^-at measures the nucleon mass shift away from the 
chiral limit and is particularly suited to test our understanding of the mechanism of 
the spontaneous and explicit chiral symmetry breaking in QCD due to the non-zero 
u, d quark masses ( For an elegant review of the earlier works, one can consult 
Ref. P ). The precise knowledge of the values of the T, nN is of great importance for 
many phenomenological applications, for example, the E^jv enters the counting rates 
in searching for the Higgs boson [2], supersymmetric particles [3] and dark matter 
|U |3] . However, no experimental method can be used to measure the T, wN directly. 
The low energy theorem relates the nucleon scalar form-factor a(t) to the isospin- 
even ttN scattering amplitude D + (u,t) at the un-physical Cheng-Dashen point, v = 
0,t = 2m 2 [Oj. The Cheng-Dashen point lies outside the physical ttN scattering 
region, we have to extrapolate the experimental D + amplitude to obtain the E„-at 
with the general techniques of the dispersion relation and partial-waves analysis, 
the bar over D + indicates that the pseudo-vector Born term has been subtracted. 
Earlier analysis performed by Koch [7] and Gasser, Leutwyler, Sainio [S] gave the 
canonical value for the a(2m 2 ), a(2m 2 ) « 60MeV, however, the recent analysis 
of the tiN scattering data supports the values E^ = 79 ± 7 MeV Although 
there have been a lot of works on the pion-nucleon sigma-term, for example, chiral 
perturbation theory [TQ1 HI1 H2] , lattice QCD [T3J EH 113 EH] , various chiral quark 
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models [T71 HH1 EH 1201, or Schwinger-Dyson equation [21], the value of the sigma 
term remains a puzzle. 

In this article, we calculate the scalar form-factor a(t) of the proton in the 
framework of the light-cone sum rules (LCSR) [221 EHI which combine the stan- 
dard techniques of the QCD sum rules with the conventional parton distribution 
amplitudes describing the hard exclusive processes [2H1- In the LCSR approach, 
the short- distance operator product expansion with the vacuum condensates of in- 
creasing dimensions is replaced by the light-cone expansion with the distribution 
amplitudes (which correspond to the sum of an infinite series of operators with the 
same twist) of increasing twists to parameterize the non-perturbative QCD vacuum, 
while the contributions from the hard re-scattering can be correctly incorporated 
as the O(ot s ) corrections |24| . In recent years, there have been a lot of applica- 
tions of the LCSR to the mesons, for example, the form-factors, strong coupling 
constants and hadronic matrix elements [213, the applications to the baryons are 
cumbersome and only the nucleon electromagnetic form-factors [211 an d the weak 
decay A& — > pivi [21] are studied, the higher twists distribution amplitudes for the 
baryons were not available until recently [2T?j . 

The article is arranged as follows: we derive the light-cone sum rules for the 
scalar form-factor a(t) of the proton in section II; in section III, numerical results 
and discussion; section VI is reserved for conclusion. 



2 Light-cone sum rules for the scalar form-factor 

In the following, we write down the two-point correlation function IT(P, q) in the 
framework of the LCSR approach, 

U{P,q)=i J d*xe iq - x {0\T{r}{0)J{x)}\P), (1) 

with the scalar current 

J(x) = u(x)u(x) + d(x)d(x), (2) 
and the baryon current [28J 

77(0) = e«* [^(0)67^(0)] 75^(0), 
(0\V{0)\P) = f V {P-z)jN{P), (3) 

here z is a light-cone vector, z 2 = 0, and the Jn is the coupling constant for the lead- 
ing twist distribution amplitude [3Uj. At the large Euclidean momenta P' 2 = (P—q) 2 
and q 2 = —Q 2 , the correlation function II(P, q) can be calculated in perturbation 
theory. In calculation, we need the following light-cone expanded quark propagator 
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i 



+ ^/L / dv \ (1 - v) ^V^H + ^/^"M jtf + ...,(4) 



where = g s G° Ml/ (X a /2) is the gluon field strength tensor and d is the space- 
time dimension. The contributions proportional to the G^ u can give rise to four- 
particle (and five-particle) nucleon distribution amplitudes with a gluon or quark- 
antiquark pair in addition to the three valence quarks, their corrections are usually 
not expected to play any significant roles [32] and neglected here [211 123 • Employ 
the light-cone quark propagator in the correlation function H(P, q), we obtain 

/ iq-x 
d'x ^ {2(C7^( 75 We« k (0\T «(0)^(x)<(0)} \P) 

+ (C/)^( 7 5/^)V^(0|T{<(0)^(0)<(x)} \P)} . (5) 

In the light-cone limit x 2 — > 0, the remaining three-quark operator sandwiched 
between the proton state and the vacuum can be written in terms of the nucleon dis- 
tribution amplitudes [HO] [2E1 [2H] • It is obviously that if we only take into account the 
three valence quark component of the distribution amplitudes, the correlation func- 
tion % J d^x e iq ' x (0\T {i](0)s(x)s(x)} \P) must be zero, i.e. the strange-component 
of the scalar form-factor (P'\s(0)s(Q)\P) = 0; if the strange-component of the itN 
sigma term manifests itself, the distribution amplitudes with additional valence glu- 
ons and quark- ant iquark pairs must play significant roles. The three valence quark 
components of the nucleon distribution amplitudes are defined by the matrix ele- 
ment, 



4(0|e ijfc <(ai^)^(«2x)<(a 3 a;)|P) = S 1 MC 0lP { l5 N) 1 + S 2 M 2 C a ^ l5 N) 1 

r 2 M 2 

+ 7W(7 5 CViV 7 + V 2 M 2 ( l5 C) a ^N), + (Vi + ^ r V 1 M )( y PC) a/3 ( 75 iV) 7 

+ V 2 M{PC) a(3 {^N) 1 + V 3 M( ltl C) a p(Yl 5 N), + V 4 M 2 (^C)^( 75 iV) 7 
+ V 5 M 2 ( lfl C) aP (ia^x ul5 N\ + V 6 M 3 ^C) a/3 (^ 5 N)^ 
x 2 M 2 

+ (Ai + -^A^){p lh C) a pN^ + A 2 M(p l5 C) a ^N), + A 3 M( W C) a p(^N^ 
+ AM 2 (^ 75 C) Q/3 iV 7 + A 5 M 2 ( lfll5 C) a p(ta^x u N), + A 6 M 3 (^ 5 C) a ^N), 

T 2 M 2 

+ (7! + ^ r -T 1 M )( J P^^C) a / 3 (7 M 7 5 iV) 7 + T 2 M(x»P^a, u CU( l5 N), 



+ r 3 M(^C) Q/3 (a^ 7 5iV) 7 + r 4 M(PV M ^) Q/3 (a w x p75 iV) 7 

+ r 5 M 2 (x^a^C) Q/3 (r 75 iV) 7 + T 6 M 2 (^P^^C) Q/3 (^ 75 iV) 7 

+ T 7 MV^) Q/3 (a^ 75 iV) 7 + T 8 M 3 (xV^CVK% 75 iV) 7 . (6) 
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The calligraphic distribution amplitudes do not have definite twist and can be related 
to the ones with definite twist as 

Si — Si, 2P ■ XS2 — Si — S2, 
Vi = Pi, 2P-xV 2 = P 1 -P 2 

for the scalar and pseudo-scalar distribution amplitudes, 

Vi = Vi, 2P ■ xV 2 = Vi-V 2 - V 3 , 

2V 3 = V 3 , 4P • xV 4 = -2Vi + V 3 + V A + 2V 5 , 

4P ■ xV 5 = V 4 - V 3 , (2P ■ x) 2 V 6 = -Vi + V 2 + V 3 + V 4 + V 5 -V 6 

for the vector distribution amplitudes, 

Ai = A u 2P ■ xA 2 = -A 1 + A 2 - A 3 , 

2A 3 = A 3 , AP ■ xAa = -2Ai - A 3 - A 4 + 2A 5 , 

4P • xA 5 = A 3 - A 4 , (2P ■ x) 2 A 6 = Ai - A 2 + A 3 + A 4 - A 5 + A 6 

for the axial vector distribution amplitudes, and 

T 1 = T 1 , 2P ■ xT 2 = Ti + T 2 — 2T 3 , 

2T 3 = TV, 2P • x% = T 1 -T 2 - 2T 7 , 

2P ■ x% = -Ti + T 5 + 2T 8 , (2P • xf% = 2T 2 - 2T 3 - 2T 4 + 2T 5 + 2T 7 + 2T 8 , 

4P • xT 7 =T 7 - T 8 , (2P • x) 2 T 8 = -T x + T 2 + T 5 - T 6 + 2T 7 + 2T 8 

for the tensor distribution amplitudes. The distribution amplitudes F = Vi, Ai, Ti, 
Si, Pi can be represented as 

F(a iP ■ x) = J Vxe- ip - x ^ x ' a 'F( Xi ) , (7) 

with 

Vx = dxidx 2 dx 3 S(xi + x 2 + x 3 — 1). 

Those light-cone distribution amplitudes are scale dependent and can be expanded 
with the conformal operators, to the next-to-leading conformal spin accuracy, we 
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The Vi, Ai and Ti are leading twist-3 distribution amplitudes; the Si, P\, V 2 , V 3 , 
A 2 , A 3 , T 2 , T 3 and T 7 are twist-4 distribution amplitudes; the 5*2, P 2 , V4, V 5 , At, 
A 5 , T 4 , T 5 and T 8 are twist-5 distribution amplitudes; while the twist-6 distribution 
amplitudes are the Vq, Aq and T 6 . Those parameters 4> 3 , 0°, 4> Q A) 0°, £4, £5, ^4, 1/;®, 

03 > 03 > 04 > 04 > ^4 , ^4~, £4 > £4", 05 > 05 > V's » ^5~> ^5~) 65") 06 > 06 can be expressed 
in terms of eight independent matrix elements of the local operators, for the details, 
one can consult Ref . . 

Taking into account the three valence quark light-cone distribution amplitudes 
up to twist-6 and performing the integration over the x in the coordinate space, 
finally we obtain the following results, 

U(P,q) = /P - zN(P) |m J dt 2 t 2 dti 1 



{q-t 2 Pf 

{2 [Si + T 7 ] (ti, t 2 , 1 - h - t 2 ) - V 3 (h, 1-ti- t 2 , t 2 )} 

pi pX rl — 12 

+ Ml dX dt 2 / dti 



,<■ ,1 Jo L (q-^P? 
{[Vi -V 2 - V 3 ] (t u l-ti- t 2 , t 2 ) - 2 [Tx - T 3 - T 7 ] (t a , t 2 , 1 - h - t 2 )} 

+ AM [ dXX [ dt 2 f 2 dh (g . ~ A f I'J [T 2 -T 3 + T 7 ] (t u t 2 , 1-h-t 

JO Jl JO [Q — XP) 



+ M 3 / drr dX dt 2 / dti 



o 



i Ji Jo (q-rPy 
[-Vi + V 2 + V 3 + V 4 + V 5 - V 6 ] (h, t 2 , t 2 )} 

+ .... (9) 

According to the basic assumption of current-hadron duality in the QCD sum 
rules approach j2S] , we insert a complete series of intermediate states satisfying the 
unitarity principle with the same quantum numbers as the current operator 7/(0) 
into the correlation function in Eq.(l) to obtain the hadronic representation. After 
isolating the pole terms of the lowest proton state, we obtain the following result, 

tP> ■ zf N N(P - q)(N(P - q)\u(0H0) + d(0)d{0)\N(P)) 

n(P ' 9) = M 2 — (q — P) 2 + " 

z*P'-zf N {P-4+M}a(t)N(P) + _ # 



in [M 2 -(q- P)' 



here m = m "+ md . The structure jt has an odd number of 7-matrix and conserves 
chirality, the structures /P , jq have even number of 7-matrixes and violate 
chirality. In the original QCD sum rules analysis of the nucleon magnetic moments 
|33j . the interval of dimensions (of the condensates) for the odd structure is larger 
than the interval of dimensions for the even structures, one may expect a better 
accuracy of the results obtained from the sum rules with the odd structure. In this 
article, we choose the structure / for analysis. 
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The Borel transformation and the continuum states subtraction can be performed 
by using the following substitution rules, 

p(x) f 1 dx p(x) f 1 dx 



/p(x) f dx p(x) f dx , . 

dx -. ' = - / 9 - / — p(x)e 
(q - xP) 2 J Q x s - P' 2 J X0 x I K ' 



p(x) f dx p(x) 1 f dx . . -t^t p(xo)e m b 



(q-xPy ' J x 2 (s-P' 2 ) 2 M 2 B J X0 x 2 ' KJ Q 2 + x 2 M 2 ' 



s = (l- x )M 2 + tL^-Q 2 , 

x 



^(Q 2 + s - M 2 ) 2 + 4M 2 Q 2 - (Q 2 + s - M 2 ) 
Xo = 2M 2 • 

Finally we obtain the sum rule for the scalar form-factor a(t = —Q 2 ), 

M 2 

a(t)f N e ^ 
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m 



f 1 d\ f x , f 1 '* 2 , 
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Jxn A J\ JO 



/x ■' 1 JO k ^-^B 
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2- r^/v rv g2+A2M2 



o ^ 2 A Jo 



/ A(1-A)M 2 + (1-A)Q 2 | _ T1 , 
exp | ^2 1 [T 2 - T 3 + T 7 ] (*i,* 2 , 1 — *i — *2j 

- 2mj^dt 2 J a d*iexp|--^-|[T 2 -T 3 + 7V](ti,t 2 ,l-*i-* 2 ) 
ml 2 /^dr r „ /" A , Z" 1 "* 2 , f r(l-r)M 2 + (l-r)Q 2 

[-V, + V 2 + V 3 + V 4 + V 5 - V 6 ] (h, 1 - h - t 2 , t 2 ) 
+ ^9° " — 977? / dX [ dt 2 [ dt\ exp 



M 2 + x 2 M 2 J 1 J 1 Z J 1 " 1 M 2 
[-V-, + V 2 + V 3 + V 4 + V 5 - V 6 ] (t ± , 1 - h - t 2 , t 2 ). (12) 



3 Numerical results and discussions 

The input parameters have to be specified before the numerical analysis. We choose 
the suitable values for the Borel parameter Mb, M b = (1.5— 2.5)GeV 2 . In this range, 
the Borel parameter Mb is small enough to warrant the higher mass resonances and 
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o Q 2 =4GeV 2 ; 
Q 2 =5GeV 2 . 
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Figure 1: The a(t) with the Borel parameter M\ for s = 2.25GeV^ 2 . 

continuum states are sufficiently suppressed, on the other hand, it is large enough 
to warrant the convergence of the light-cone expansion with increasing twists in 
the perturbative QCD calculations [SHI E3| The numerical results show that in 
this range, the scalar form-factor a(t = —Q 2 ) is almost independent on the Borel 
parameter Mg, which we can see from the Fig.l for Q 2 = 3GeV 2 , AGeV 2 and 5GeV 2 . 
In this article, we take the middle point for the Borel parameter, M\ = 2.0GeV 2 in 
numerical analysis, such a specialization will not lead to much uncertainties on the 
final results and impair the predictive ability. 

There are two independent interpolating currents with spin ~ and isospin ~, both 
are expected to excite the ground state proton from the vacuum, the general form 
of the proton current can be written as [37] 

J(x,t) = e abc {[ul(x)C l5 d b (x)]u c (x)+t [ul(x)Cd b (x)] l5 u c (x)}, (13) 

in the limit t = —1, we recover the Ioffe current. The Monte- Carlo calculations 
for the two-point vacuum correlation function indicate that the optimal mixing be 
t = —1.2, the threshold parameter be ^/sq = (1.53 ± 0Al)GeV and the mass of the 
proton be M = (1.17 ± 0.26)GeV jSH], furthermore, the Monte-Carlo method has 
been successfully applied in studying the axial coupling constant and the magnetic 
moments of the decuplet baryons j^Sj. Here we take the Ioffe- type current r)(x) in 
Eq. (3) to keep in consistent with the sum rules used in determining the parameters 
in the light-cone distribution amplitudes, ^3, 0g, 0°, 0°, £4, £5, ipl, ip®, 

07 , (j>t, ip±, ipti 6T » £l> ^5 > <t>t , ^5", ^t, is, \ti ^6 > ^6 I furthermore, we take 
the physical mass for the proton rather than determine it from the corresponding 
sum rules, this obviously leads to some deviations from the Monte-Carlo adjusted 
threshold parameter s . In the Fig. 2, we plot the dependence on the threshold 
parameter s f° r the sum rules in a large range, ^/sq = (1.5 ± 0A)GeV. From 
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Figure 2: The a{t) with the central values of the input parameters. 



9 







4>7 


4>t 






^ + 


£9 


5" 




twist-3 


i = 3 


0.53 


2.11 


0.57 














twist-4 


i = 4 


-1.08 


3.22 


2.12 


1.61 


-6.13 


0.99 


0.85 


2.79 


0.56 


twist-5 


% = 5 


-1.08 


-2.01 


1.42 


1.61 


-0.98 


-0.99 


0.85 


-0.95 


0.46 


twist-6 


% = 6 


0.53 


3.09 


-0.25 















Table 1: Numerical values for the parameters, the values are given in units of 
10" 2 GeV 2 HI]. 

the figure, we can see that in the region Q 2 > 3GeV 2 , the scalar form-factor is 
insensitive to threshold parameter sq, while in the region Q 2 < 3GeV 2 , the curves 
for the form-factor a(t = —Q 2 ) with ^/sq < 1.5GeV are quite different from those 
ones with y/so > 1.5GeV, which may be due to the contributions from the high 
resonances and continuum states. For simplicity, we choose the standard value for 
the threshold parameter sq, Sq = 2.25GeV 2 to subtract the contributions from the 
higher resonances and continuum states, i.e. we restrict the range of integral to the 
energy region below the Roper resonance (iV(1440)); on the other hand, it is large 
enough to include all the contributions from the proton. The current masses for the 
u and d quarks are, m u = (1.5 — A)MeV and = (4 — 8)MeV from the Particle 
Data Group in 2004 [33], we choose m = = 5MeV. For Q 2 = (2 - 5)GW 2 , 

x > Xq = 0.5 — 0.7, and the average value (x) = 0.75 — 0.85, with the intermediate 
and large space-like momentum Q 2 , the end-point contributions (or the Feynman 
mechanism) are dominant 2 , it is consistent with the growing consensus that the 
onset of the perturbative QCD region in exclusive processes is postponed to very 
large energy scales. 

The parameters in the light-cone distribution amplitudes 0°, 0g, 0°, 0°, 

1>l, ^5> 03> 03"> <f>4, <f>t, 1>Z, V>4 > €i, 05 > <f>t, ^5, ^t, £t , <f>6 , 06 are 

scale dependent and can be calculated with the corresponding QCD sum rules, the 
approximated central values are presented in the Table 1. They are functions of eight 
independent parameters, f^, Xi, A 2 , Vf, A™, /J, /J and f*, the explicit expressions 
are presented in the appendix, for detailed and systematic studies about this subject, 
one can consult Ref . j2Hl • Here we neglect the scale dependence and take the following 
values for the eight independent parameters, /jy = (5.3 ± 0.5) x 10~ 3 GeV 2 , Ai = 
-(2.7 ± 0.9) x 10~ 2 GeV 2 , A 2 = (5.1 ± 1.9) x lO~ 2 GeV 2 , Vf = 0.23 ± 0.03, A\ = 
0.38 ± 0.15, ff = 0.6 ± 0.2, f$ = 0.15 ± 0.06, = 0.22 ± 0.15. In estimating 
those parameters with the QCD sum rules, only the first few moments are taken 
into account, the values are not very accurate. One can map the uncertainties 
of the input parameters into the uncertainties of the adjusted phenomenological 
parameters (331^3, for example, the threshold parameter s , the mass of the proton 
M, the scalar form-factor a(t = —Q 2 ), etc, with the Monte-Carlo method through 

2 Our work on the axial form- factor of the nucleons with the LCSR also leads to the same 
observation, and will be presented elsewhere. 
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Figure 3: The a(t) with the Q 2 for the central values of the parameters except Ai . 

\ 2 minimization, which may be the most realistic estimates of the uncertainties as 
there are many input parameters which can be taken into account simultaneously, 
however, we are no expert in Monte-Carlo simulation, the traditional uncertainties 
analysis is chosen in this article. 

We perform the operator product expansion in the light-cone with large Q 2 and 
P /2 , the scalar form-factor a(t = —Q 2 ) make sense at the range Q 2 > 2GeV 2 , with 
the low momentum transfers, the operator product expansion is questionable. In 
this article, we devote to calculate the scalar form-factor at the range Q 2 > 2GeV 2 , 
which corresponding to the size about 0.1 fm, (q — xP) 2 —>■ xM\ after the Borel 
transformation, in the region Mr = (1.5— 2.5)GeV 2 , J-— < _ = , „ .. ~ 0.24 fm, 

' to V ii / xM 2 — V0.5xl.5GeV J ' 

retaining only the three valence quark light-cone distribution amplitudes up to twist- 
6 is reasonable. The size of the proton is about the order of lfm which corresponding 
to the confinement scale Aqcd ~ 0.2GeV, we only investigate short distance physics 
inside the proton. With smaller momentum transfers, the contributions from the soft 
(small virtual) gluons and quarks become larger, and the multi-parton configurations 
become more and more important, the quark and gluons degrees of freedom have 
to be integrated out, we can work in the hadronic representation and resort to the 
chiral perturbation theory to deal with the problems |10l ITT| IT2"] . In numerical 
analysis, we observe that the scalar form-factor a(t = —Q 2 ) is sensitive to the four 
parameters, Ai, ff, /| and which are shown in Fig.3, Fig.4, Fig. 5 and Fig.6, 
respectively. Small variations of those parameters can lead to large changes 
for the values, the large uncertainties can impair the predictive ability of the sum 
rules, and those parameters Ai, ff, f$ and should be refined to make robust 
predications. The final numerical values for the scalar form-factor a(t = —Q 2 ) at 
the intermediate and large space-like momentum regions Q 2 > 2GeV 2 are plotted 
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Figure 4: The a(t) with the Q 2 for the central values of the parameters except ff. 
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Figure 5: The a(t) with the Q 2 for the central values of the parameters except f£. 
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Figure 6: The a(t) with the Q 2 for the central values of the parameters except 

in the Fig. 7, from the figure, we can see that the values of the scalar form-factor 
a(t = —Q 2 ) are compatible with the calculations of lattice QCD ^3] and chiral quark 
models [36J. If we take the values from the recent analysis of the ttN scattering data, 
er(0) pa = 79 ± 7MeV), as input, the results from the lattice calculation give 

a(-2GeV 2 ) < <r(0) x 0.1 pa 7.9MeV, our results \a(-2GeV 2 )\ pa (2.9 ± 2.7)MeV 
are smaller, however, reasonable. The a s corrections to the scalar form-factor of 
the proton may be significant, quantitative conclusion can be reached after the 
solid calculations, the calculations are tedious though not impossible, and beyond 
the present work. In the case of the ir meson, the a s corrections of the twist-2 
light-cone distribution amplitude reproduce the behavior for the electro-magnetic 
form-factor with large Q 2 , which corresponding to the hard re-scattering mechanism 
[2"lj . The contributions from the four-particle (and five-particle) nucleon distribution 
amplitudes with a gluon or quark-antiquark pair in addition to the three valence 
quarks are usually not expected to play any significant roles [32] and neglected here. 
The consistent and complete LCSR analysis should include the contributions from 
the perturbative a s corrections, the distribution amplitudes with additional valence 
gluons and quark-antiquark pairs, and improve the parameters which enter in the 
LCSRs. 

4 Conclusion 

In this work, we calculate the scalar form-factor a(t = —Q 2 ) of the proton in the 
framework of the LCSR approach up to twist-6 three valence quark light-cone distri- 
bution amplitudes and observe the scalar form-factor a(t = —Q 2 ) with intermediate 
and large momentum transfers, Q 2 > 2GeV 2 , has significant contributions from 
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Figure 7: The a(t) with the Q 2 . 



the end-point (or soft) terms, it is consistent with the growing consensus that the 
onset of the perturbative QCD region in exclusive processes is postponed to very 
large energy scales. In numerical analysis, we observe that the scalar form-factor 
a(t = —Q 2 ) is sensitive to the four parameters, Ai, ff, f$ and small variations 
of those parameters can lead to large changes for the values. The large uncertain- 
ties can impair the predictive ability of the sum rules, the parameters Ai, ff, f$ 
and fi should be refined to make robust predications. The numerical values for 
the a(t = —Q 2 ) are compatible with the calculations from the chiral quark model 
and lattice QCD. The consistent and complete LCSR analysis should include the 
contributions from the perturbative a s corrections, the distribution amplitudes with 
additional valence gluons and quark-antiquark pairs, and improve the parameters 
which enter in the LCSRs. 
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Appendix 



€2 = ^ = |a 2 , ^ = ^ = \{f N -\i) ■ 

03 = 

3 + = ^(l-3Vf), 

04" = ^ (Ai(l - 2ft - 4/r) + f N (2A u 1 - 1)) , 

4 + = | (Ai(3 - 10/i) - /^(lOV^ - 3)) , 

C = (Ai(2 - 7/? + + f N (A? + W? - 2)) , 

= -\ (Ai(-2 + 5/? + 5/f ) + M2 + 5A? - 5Vf)) , 

C = ^A 2 (4-15/ 2 d ), 

= ^A 2 (4-15/ 2 d ), 

05 = |(Ai(y?-/r) + M2A?-l)) , 
5 + = -^(A 1 (4/ 1 d -l) + / JV (3 + 4y i d )) , 

C = ^(Ai(A d -/r) + M2-^-3^)) , 

= ~\ (Ai(-1 + 2/f + 2tf) + /at(5 + 2A« - 2V?)) , 

£5 = _ ^A 2 / 2 i ) 

e 5 + = ^A 2 (2-9/ 2 d ), 

06 = 5 (Ai(l - 4/* - 2tf) + Ml + 4A?)) , 
6 + = -^(1- 2/^ + ^(4^-1)) . 
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